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1 Introduction 

Since the pioneering works pQ-[2] the Toda hierarchy has become one of the paradigmatic examples 
of the relevance of integrable systems in the theory of random matrix models. As a consequence 
of the activity in this field a rich theory of the different facets of the Toda hierarchy has been 
developed. 

The present work is motivated by the applications of the Toda hierarchy theory to the Hermitian 
matrix model. In this model the first integrable structure which emerges is the 1-Toda hierarchy 



on semi-infinite tridiagonal matrices 

L = A + u n + v n A T , n>0. 

Here A is the standard shift matrix and A + denotes the upper part (above the main diagonal) of 
semi-infinite matrices A. This relationship may be conveniently described by considering infinite- 
dimensional deformations of monic orthogonal polynomials on the real line 

P n {z,t) = z n + --- , t:={tx,t 2 ,...), n>0, 

with respect to an exponential weight: 

/oo 
P n (z,t)P m (z,t)e v ^ t Uz = h n (t)5 nm , V(z,t) :=^2(t k + c k ) 
-°° k>l 



z\ 



where c := (ci, C2, . . .) is a given set of complex constants. It turns out (see for instance [3]) that 
the functions 

ip n (z,t) := P n (z,t) exp(^ z k t k ), n>0, 

k>l 

satisfy the linear system of the semi-infinite 1-Toda hierarchy 

^- = (L k ) + Tp n , Lti) n = ztp n , k>l, n > 0, 
dt k 

and have a r-function representation 

provided by the partition function of the Hermitian matrix model 

N 



TN {t) = Z N (t):= f f[(dx k e v ^yA( Xl ,--- ,x N )) 2 , (1) 
J RN k=i 

where A(xi, ■■■ ,x n ) := Ili^^t _ x j) 
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Many exciting properties of the Hermitian model emerge in the analysis of its large iV-limit 
j3)-[H]- One of the main tools supplied by the theory of the Toda hierarchy for such analysis PQ is 
the use of a pair of constraints called string equations 

L = A + u n +v n A T , M = ^2k{t k + c k )(L k ~ 1 ) + (2) 

k>l 

which are satisfied by the canonically conjugated operators 

Ltp n = z ip n , = M ip n . 

oz 

The present paper deals with the analysis of the large iV-limit of the partition function 

f N 

Z N (Nt)= / H(dx k e NV ^))(A( Xl ,--- ,x N )) 2 . (3) 
k=i 

Here a small parameter e := 1/N and rescaled variables t := et and constants c := ec have been 
introduced. Since e plays the role of the Planck constant h, expansions in powers of e are referred 
to as semiclassical expansions. The same slow variables t = e t together with a continuous variable 
x := en are introduced to pass from the Toda hierarchy to its dispersionful formulation [10], which 
provides an interpolated continuous version of the Toda hierarchy. In this way, and due to the fact 
that T n (t) = Z n (t) is a r- function of the semi-infinite 1-Toda hierachy, it is natural to expect that 
a r-function r(e,x,t) of the dispersionful 1-Toda hierarchy verifying 

T(e,en,t) = Z n (Nt), (4) 

should describe the large iV-limit of the Hermitian model. In this paper we are concerned with the 
characterization of this solution of the dispersionful 1-Toda hierarchy. 

The main result of our work is a scheme for obtaining solutions of the dispersionful 2- Toda 
hierarchy satisfying the system of string equations 

£ = £, M + F(£)=M+F(£), (5) 

where (£,A4) and (£, Ai) denote two pairs of Lax-Orlov operators and (F, F) are two arbitrary 
functions. The first string equation represents the 1-Toda reduction condition and is satisfied by 
Lax operators of the form 

C = £ = A + u + uA" 1 , 



where now A := exp (ed x ), and (u, v) are characterized by semiclassical expansions 

J2e 2h v m - (6) 



fc>0 fe>0 



The point is that for x = 1 and F = the constraints © interpolate ([2]), so that the corresponding 
solution of the dispersionful 1-Toda hierarchy is a candidate to the solution underlying the large 
iV-limit of the Hermitian model. That it is the only possible candidate can be argued as follows: 



1. Recent research [12j-|14j proved that the solutions of an extended version of the dispersionful 
1-Toda hierarchy are determined by the leading order terms (u^°\ v^). In fact, the coeffi- 
cients (u^ k \ t/ 2fc J) are rational functions of (u^ \v^ ') and their x-derivatives (quasi-triviality 
property). 
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2. As it is shown in this paper, the terms (u^°\v^) of the solution of ([5]) coincide with those 
characterizing the leading order in the large A^-expansion (planar limit) of the Hermitian 
model. 

Our strategy is inspired by previous results |15j-|18j on solution methods for dispersionless 
string equations. We also develop some useful standard technology of the theory of Lax equations 
|19|-|21j in the context of the dispersionful 1-Toda hierarchy. Thus we introduce two generating 
functions M and T related to the resolvent of the Lax operator which play a crucial role in our 
analysis. 

The paper is organized as follows: 

In the next section the basic theory of the dispersionful 2-Toda hierarchy and the method of 
string equations are discussed. In Section 3 we deal with the dispersionful 1-Toda hierarchy and 
its relationship with the Hermitian matrix model from the point of view of the continuous string 
equations ©. The generating functions R and T are introduced and are characterized by two 
important identities. Our main results are derived in Section 4 where a scheme for solving the 
string equations © in terms of semiclassical expansions is provided. In particular we prove that 
the leading terms of these expansions characterize the planar limit of the Hermitian matrix model. 
In Section 5 it is showed how the double scaling limit method can be naturally implemented in our 
scheme. 

Applications of our method to normal matrix models which are also related to string equations 
of the Toda hierarchy [22]- [29] will be considered elsewhere. 

2 String equations in the dispersionful 2-Toda hierar- 
chy 

2.1 The dispersionful 2-Toda hierarchy 

The formulation of the dispersionful 2-Toda hierarchy [10] uses operators of the form 

A = '^2a j (e,x,t,i)A j , A := exp (e d x ) , (7) 

where x is a complex variable and the coefficients are in turn series in the small parameter e 

aj (e,x,t,t) = ^e fc af } (x,t,t). 

fcez 

Here t = (ti, t%, . . .) and t = (ti, tjj) • • •) denote two infinite sets of complex variables. The order in 
e of A is defined by 

ord e (^) := max{-/c| af\x,t,i) / 0}. 

For example ord e (e) = —1 and ord e (A) = 0. In particular, zero-order operators are those with 
regular coemcents a,j as e — > 0. As usual the A± parts of A will denote the truncations of A-series 
in the positive and strictly negative power terms, respectively. Given a function w depending on 
x, the following notation convention will be henceforth used 

u)[ r ] := A r w = w(x + r e), r £ Z. 
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The dispersionful 2-Toda hierarchy can be formulated in terms of a pair of formal wave functions 
of the form 



1 °° 1 

* = exp (~ s), § = h zj +x\ogz~Y^ —Sj+i, 



7 

3=1 3>1 



(8) 



1 00 1 

^ = z~ l exp ^— §J , S = ^2 tj^ ~ x 1°§ z ~ $o — ^2 ~ T ~~jSj+i, 

6 3=1 3>i 3 Z 

where 

s i = E e " 5 f ( x ' *' *)> ^ = E efc ( x > *> *)> 

fc>0 fc>0 

These functions = ty, are assumed to satisfy the linear system 

where the Lax operators £ and £ are determined by the equations 

£^ = zt>, £~$ = z^, (10) 

and are assumed |10j to be of zero order in e. We will also use the Orlov operators Ai and A4 
characterized by 

d^f dW 

M* = e—; M* = e—, (11) 

which satisfy 

[£,M] = [£,M] =e. 
Using ([S]) and (|1U|) - (|11|) one sees that the following expansions follow 

00 

£ = A + n + ui A' 1 + • • • , M = 3 tj£P~ X + x£- x + Y^ S j+1 £~i~ x 

3=1 3>1 

(12) 

00 

£ = u_i A" 1 + u + «i A + • • • , M = Y J - + e)£- 1 + Y S3+1 

3=1 3>1 

Furthermore, Q can be rewritten in Lax form as 

t f r U')„K], .™-l<P)-,K\, (13) 

where K = £, M, £, M. 



5 



There is also a r-function representation of the wave functions [TO] 

/l ,v~^ i , x\ r(e,x,t — e [z _1 l,t) 
=exp(- YYV+xlogz ) [ rT jLj ; 



i=i 



(14) 



where [2 x ] := (l/z, l/2z 2 , l/3z 3 , . . .) and r is of the form. 

r = exp(iF), F = ^e fc FW(x,t,t). (15) 
' fc>0 

The dispersionful 2-Toda hierachy arises as a continuum limit of the standard 2-Toda hierachy 
[9] in which the standard discrete variable n is substituted by a continuous variable x and two sets 
of fast continuous variables t := e _1 t, i := e _1 t are introduced. Thus, the main dynamical objects 
(r-functions, wave functions, Lax and Orlov operators) of both hierarchies are related by 

T(e,en,et,et) = T n (t,i) (16) 

and 

e, en, et, et) = t/j n (z,t,t), *f>(z,e,en,et,ei) = $ n (z,t,i), 
£(z,e,en,et,ei) = L(z,n,t,i), C{z,en,et,ei) = L(z,n,t,i), (17) 
Ai(z, e, en, et, ei) = e M{z, n,t,i), M.(z,e,en,et,ei) = e M(z, n, t, i). 

Our subsequent analysis uses an important result proved by Takasaki and Takebe (Proposition 
2.7.11. in [TO]) 

Theorem 1. Suppose that 

V(e,xA-\A) = ^ftfcxA- 1 )^, Q(e,xA" 1 ,A) = ^(e^A" 1 ) A fc , 

V(e,xA-\A) = Y,Pk(e,xA- 1 )A k , Q(e,xA-\ A) = J^ftfoaA" 1 ) A k , 

are operators of zero order in e verifying 

[V,Q] = [V,Q]=e 

If (C,M.) and (C,M.) are operators of zero order in e of the form ([TO]) which satisfy the pair of 
constraints 

V(e,M,£)=V{e,M,C), Q(e,M,£) = Q{e,M,C), (18) 

and 

[£,M] = \C,JZ\ = e, (19) 
then (jO,A4) and {C,M) are solutions of the Lax equations (I13p of the 2-Toda hierarchy. 
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Constraints of the form (|18p are called string equations. In this work we are interested in the 
particular example given by 

= (20) 

[M + F{£) =M + F{£), 

where F{C) and F(C) are arbitrary functions of the form 

3 The dispersionful 1-Toda hierarchy and the Hermi- 
tian model 

The first string equation in ([20]) represents the so called tridiagonal (1-Toda) reduction of the 
dispersionful 2-Toda hierarchy and implies the following form of the Lax operators 

C = I = A + u + vA- 1 . (21) 

Thus, as a consequence of the Lax equations, u and v depend on (t,t) through the combination 
t — t. Moreover (|2ip implies 



(A + u + vA' 1 )^ = zV, {A + u + vA~ 1 )^ = z^, (22) 

so that 

u = e- 1 (S 2[1] -S 2 ), logv = e- 1 (S ohl] -S ). (23) 



In order to solve the string equations (|20p it is required to characterize the action of the 
operators (C J )+ and (C 3 )- on the wave functions *Sf and ^. This calculation is also needed to 
determine the integrable systems of the dispersionful 1-Toda hierarchy. We start by introducing 
the two series in z 

p(z) = z-u + ofy, p( z ) = ^ + 0^, z^oo, (24) 

satisfying 

A$>=p(z)fy, A* =p(z)"3>, (25) 
which according to (|22p are determined by 

p(z)+u + Vt = z ' ( 26 ) 

P[-i]{ z ) 

where p = p, p. By using ([22]) it is clear that there are functions ay, Pj,otj, j3j , which depend 
polynomially in z, such that 

e — = (tf )+¥ = aj * + p j A* = { aj + f3 3 p)*, 

(27) 

e ^ = = a, * + 0j A* = (a, + ft p)*, 



and 



Hence, we have 



j — z° ccj, f3j — Pj. 



otj + Pj P = a.s(z) = z j + o 



a 3 +P j p = d t Mz) = -d tj So + 



so that 



1 



a j = n \ z ° ~ d h S ® ~ 0^' (P + P))® ' #7 



p — p 



(28) 
(29) 
(30) 



where ( )® and ( ) e stand for the projections of z-series on the subspaces generated by the positive 
and strictly negative powers , respectively. 

At this point it is useful to introduce the generating functions 



E 



Rk(u,v) 



T := 



p+p 



p — p f — ' z" 

1 * fc>0 



E 



T k (u,v) 



p — p L — ' z" 



Ro — To — 1. 



By substituting p and p by their expressions in terms of R and T in the identities 

p-p 



PP[-i] ~PP[-i] 
u = z H ^ — - - — -, v 



P[-l] ~ P[-i] 



P[-i) ~ P[-i] 



P[-l]P[-l], 



(31) 



(32) 



we obtain the following relations 



T [1] +T+-(u [l] -z)R [l] =0, 



(33) 



T 2 



I"[i]«[i] =1 . 



which allow us to compute recursively the coefficients of the series (|3ip as polynomials in u, v and 
their x-translations u\ r ] and ur r i . Indeed, the system (l33j) implies 



2T k+ i = - ^ Ti Tj + 4u[i] ^ 

^fc+i = ^-Rfc + gPfc+i + r fc+i[-i]]' • 

For example, the first few coefficients are: 

Ti =0, i?i = u, T 2 = 2v[ X ] , R 2 = u 2 + V[x] + v, 

T 3 = 2v[ 1 j(u + U[!]), i?3 = u 3 + 2it«[ 1 ] + 2-ut; + «[i]«[i] + 

T 4 = 2«[!] (U[!j + + U 2 + U[ 2 ] + + u). 

In this way, by taking into account the second equation of (I34p . one finds 

d t . E(z) = a, +(3 jP = j - \d t] S - (z^ : 



(34) 



>e + (si 



(35) 



^-^(^5o + ^)-^ + i[-i] + o(^ 



8 



so that 

dtjS = -Rj, d tj S 2 = - 
and then from (|23p we get that the flows of the dispersionful 1-Toda hierarchy can be expressed as 

ed h u = 2 ( T J' +1 - T i+i[-i])' e^.v = v(i? i (36) 
Furthermore, our calculation implies the following useful relations 

1 „ z 



^ = (j* + 21 ^ R » e + (S5 T R) »)e) »■ 



(37) 



for j > 1. In particular, by taking the second equation of (|34p into account one finds that as 

z — ► 00 

(^)-*=(i r i + n-i] + (F))*. 

(38) 

1 „ „ / 1 \ \ _ 

2z 



We observe that since Ro = 1, Ti = these last equations hold for j > 0. 

By following the analysis of [21] it can be seen that K and T are closely related to the resolvent 
of the Lax operator C 

K:= (z-C)-\ 



Thus, from Lemmas 3.5 and 3.18 of [21] one proves that 

21 



, A ft+=R, Res7£+ = -, 
(1+T) J + + z 

where Res(][] A fe ) := Co 
r-function representation 

It follows from (|14p and (|22p that the functions u and v can be written in terms of the r-function 

as 

d r(e,ac + e,t) r(e, x + e, t) r(e, sc - e, t) 
U = 6 ^ lQg r( £ ,,,t) ' " = T^t) ' (39) 

where we have set t — t — > t. On the other hand, it can be proved p2]-[T3] that the e-expansion of 
the r-functions of the dispersionful 1-Toda hierarchy is of the form 

r = exp(4lF), F = Y, e2k F(2k) - ( 4 °) 
( fc>0 
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As a consequence u and v can be expanded as 

u = ^e k u^ k \ v = Y,e 2k v {2h) . (41) 

fe>0 fc>0 

Let us introduce the reduced § and M functions 

i>i J 

From (|14p we see that 

F(e, x, t - e [z -1 ]) - F(e, x, t) = e S r (e, z, x, t), (42) 
and by differentiating this equation with respect to z we obtain 

^2^-^-¥(e,x,t) = M r (e,z,x,t + e[z- 1 }). (43) 
i>i 2 3 

This identity can be rewritten in in a more convenient form. Indeed (|42p implies 

S r (e, z, x, t) — S r (e, z, x, t — e [z' -1 ]) = §y (e, z', x, t) — § r (e, z', x, t — e [z -1 ]), 
and by differentiating with respect to z and then taking the limit z' — >• z one finds 

M r (e,z,x,t - e [z" 1 ]) = M r (e, z, x, t) + e ^ j+1 r (e,z,x,t - e [z" 1 ]). 
Thus becomes 

3.1 The Hermitian matrix model 

Let us write the partition function of the Hermitian matrix model in terms of slow variables t := e t, 
where e = 1/N 



JT (dx k e NV ^)) (A(x l5 • • • , x n ))\ V(z, t) := V(t fc + c fc ) 
- : " fc=i V 7 fc >i 



z fe . (45) 



The large A^-limit of the model is determined by the asymptotic expansion of Z n (N t) for n = N 
as N — > oo 

- AT 

^iv(iVt)= / JJ^Xfce^^^^CACzi,---,^)) 2 , (46) 
•^ RJV fc=i 

It is well-known [3] that Z n (i) is a r- function of the semi-infinite 1-Toda hierachy , then in view 
of (|16p we may look for a r-function r(e, x,t) of the dispersionful 1-Toda hierarchy verifying 

r(e,en,t) = Z n (iVt) J (47) 
10 



and consequently 

r(e,l,t) = Z N (Nt). (48) 

The point is that for 

x = 1, tj = cj = 0, j > 1, (49) 

the system (f2"0"P of continuous string equations interpolates the discrete system ([2]). Hence the 
solution of the dispersionful 1-Toda hierarchy provided by (|20p can be expected to correspond 
to the r-function verifying (|47p and, as a consequence, to describe the the large iV-limit of the 
Hermitian matrix model. 

We may express the 1 /^-expansions of the main objects of the hermitian matrix model in terms 
of objects in the dispersionful 1-Toda hierarchy. For instance, from (|44p the one-loop correlator 

wiz) =-y J-<^> =-+—Y— 9 J^m 

w N zj +i \" 1V1 I z + N 2 2^ zj +i dt 

j>0 j>l 

becomes 

W (z) = - + M r (e, z, 1, t) - e £ z J+i~Qf~( e ' z ' ( 50 ) 

Loop correlators of higher order can be obtained from W(z) by application of the loop-insertion 
operator d/dV(z) [6] 

w{zi -- Zs) = d^)---dVW) W{zi) 
d V—!— 8 



dV(z) ' ZJ+ 1 dtf' 

1 ; j>i 3 



4 Semiclassical expansions 



We now turn to the solutions of the system of string equations (|20p . The first equation is solved 
by setting 

C = C = A + u + vA- 1 , 

which is in agreement with the asymptotic form (|12p required for C and C 

Let us consider the second string equation of (|20p . We look for solutions M and M verifying 
asymptotic expansions of the form (|12p . To this end we first set 



M + F(C) =M+ F{C) = Y,J (tj + c j) + E J & + ^ 

which, taking into account the first string equation, leads to 

oo oo 



(51) 



OO OO 

M = Y^jtj O- 1 - Ei ((*; + cj) - (t 3 + c,-)) (^'- 1 )+. 
i=i i=i 
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In order to satisfy (|12|) and (|19p we introduce auxiliary functions of the form 



1 °° 1 
* = eX p - ( ^2 tjz 3 + a; log z - ^ — S^+i), 



3=1 3>1 J 



(52) 



1 00 1 
q? = exp - ( ^2 ij z3 - {x + e) \ogz - Sq - —jSj+i 



and impose 



3=1 3>1 



dz 



d^f 

C^ = z^, M^ = e—. 

oz 



(53) 



Our aim is to determine u,v,A4 and M from (|53p . Now, with the help of (|38p . we have that the 
equations (f53j) for the Orlov operators read 



1 11 

i + E ij^ = E 3 (& (*i + <*)) (3- r ihl] + o( F 

(54) 

1 -1 00 11 

- + E ^ = - + ^ - + (V 1 + ^ + °(f 

i>2 j=i 

Matching the coefficients of z^ 1 in both sides of these two equations provides the same relation. 
Another relation is supplied by identifying the coefficients of the constant terms in the second 
equation of (|54p . Hence we get a system of two equations to determine (u, v) 



E£i j [(ij +c j )- (t 3 + cj)) Rj-! = 0, 

i (55) 

L 2 E^lii + q?) - (^ + c,-)) r ihl] = x. 

By equating the coefficients of the remaining powers of z in (|54p we characterize the functions Sj 
and S 1 ,- for j > 1 in terms of (u, v). Moreover, as it is proved below, the solution (u,v) provided 
by (|55p is of the form 

u = J2 e fc u {k) (x, t, t), v = J2 efc wW (*> *. *). 

fc>0 fc>0 

with t>( 2fc+1 ) = 0, Vfc > 0. Thus, by solving (f55l) we characterize operators (£,.M) and (C,M) 
which satisfy (|2U|) and all the requirements of Theorem 1. Therefore, they are solutions of the Lax 
equations for the dispersionful 2-Toda hierarchy. 

We observe that, as it is noticed by Takasaki and Takebe in [TU], solving the system of string 
equations (|20p does not determine the coefficient Sq in (|52p and therefore it does not determine a 
wave function of the dispersionful 1-Toda hierarchy. 



12 



4.1 An iterative scheme for determining (u,v) 

It is convenient to write (1551) in the form 



l-^-U z W{z) =0, f^U z T(z) = -2(x + e) 



where U denotes the function 



U (z, t, t) := (ftj + Cj) ~ (tj + cj] 
3=1 



(56) 



(57) 



and 7 is a large positively oriented closed path. Now by using the first identity of f|33|) and the 
two equations of (|56l) we find 



dz 
2~kI 



U z (T + T[_ 



dz 



TTl Z 



(z - u) U z 



dz 



-4x-2e, 



TTl 



so that (j55f) reduces to a pair of equations involving R only 

dz 



2ni z 



U z (z)R(z) =0, 



j> ^U z (z)R(z) = -2 
These equations together with the system 



(58) 



x — e. 



T [1] +T + -(u [1] -2)R [1] =0, 



T2 - ? "[i]«[i] =1 . 



(59) 



give rise an iterative scheme for characterizing (u, v) as Taylor series in e 

u = ^2e k u {k) (x,t,i), v = ^2e k v^ k \x,t,i). 

k>0 k>0 

The first step of the method is to determine the expansions 

R( z ) = e k R {k) , T{z) = ^e k , 



(60) 



fc>0 



fe>0 



in terms of (u,v). It can be done by equating the coefficients of powers of e in (|59p . Indeed, the 
coefficients of e° leads to 



((z - u(°)) 2 - 4v(°)) 3 



r (o) 



z — u 



(0) 



((z - u(°)) 2 - 4«(°)) 



2 



(61) 
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and the coefficients of e l (I > 1) yield the following system 



T^-(z-u^) 



1 



2 . v j 

i + j = i 

3 > 1 



(62) 



T (o) T (i) _ 4v (o) 



R(0) R(l) 



z z 



E (E.>^( E > 



i + j + k = l i\+i2=i 



ki + k 2 = k 
fci < Z 



1 £ T « T 



(J) 



i+j = l 
i,j > 1 

Some comments concerning these formulas are in order 

i) The equations (|62p determine each pair (TW, flW/z) from (T^, with j = 0, 1, 1. 

ii) The equations (|62p are linear with respect to T^, R^/z. Moreover, by taking into account 

(I6ip . we see that the determinant of the coefficients of and R( 1 '/z in (j62H is 

[(^_tx(°)) 2 -4i;( )]5. 

Hence it follows that the functions j z can be written as linear combinations of 

1 

r = 1,2,...,/ + ! 



((z--u(o)) 2 -4u(°)) 



((z - u(o))2 - 4u(°)) 



with coefficients depending on u^' and v^' with j = 0, 1, . . . , I and their ^-derivatives only. 

Now let us go back to the system f|58[) and find (u, v). By substituting the e expansion of R in 
([58]) we get a system of two equations for each 



dz 
2iri z 



U g (z)R®(z) = 0, 



(63) 



V^7 



dz 
27d 



U x (z)R®(z) 



-2x5, 



10 — °n, 



which determine each pair (vft\v®) recursively. Furthermore, we can eliminate the explicit de- 
pendence on (x,t,t) in the corresponding expressions since, by differentiating with respect to x 
the equations ([63]) for / = 



f 1 



1 



dz- 



7 ((* - u(°))3 - 4«[°]) 



0, 



dz- 



zU z 



(64) 



2ni ./ 7 ^ _ u (o))2 _ 4u (o)) 
14 



-2x, 



all the integrals of the form 



^ dz -, — : J) dz - j-, (65) 



27Ti -h ((s- U (0))2_ 4u (0)) r+ 2 27Ti/ 7 ^ _ u(0 ))2 _ 4^(0)^+5 

can be expressed in terms of (u^°\v^) and their x-derivatives. We observe that (j65H are the 
variables introduced in [5] to determine the large iV-expansion of the hermitian matrix model. 

Some important relations among the coefficients of the semiclassical expansions under consider- 
ation are found by realizing that given a solution (u(e, x), v(e, x), R(e, z, x), T(e, z, x)) of (f58j) - (f5"9"]) . 
then 

u(e, x) := u(— e, x + e), v(e, x) := e, x), 

R(e, z, x) := R(-e, z, x + e) f (e, z, x) := T(-e, z, x + 2e), 

satisfies (|58p -(|59 p as well. Thus, since the solution of (|58p - (|59p is uniquely determined by («(°), u(°)) 
we deduce that 

u(e,x) = u(e, x), v(e,x) = v(e,x), 



l(e,z,x) = R(e, z,x), T(e,z, x) = T(e, z,x). 



Hence we find 



1 2j ~ l ( l)^ 1 

t (2i-i) = i y iziL ^i-i-fc) ,.<^-n 

2 ^ k\ x 

fc=l 



B^- 1 ' = 0, 



(66) 



2 ^ fc! ' 2 ^ fc! 

fc=i fc=i 



for j = 1,2,.. 



4.1.1 Examples of calculations 

Using (|66p for j = 1, it is immediately found that 



r£ 3 , i 3 , (,"'J 



2 ((z - n(°)) 2 - 4i>(°)) 5 ((z - u(°)) 2 - 4t;( )) 

and 

u^ = l -u^\ v^=0. (68) 
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With the help of Mathematica, one obtains 

R (2) 4(z _ u (0) )u (2) + 8v (2) + u (0) 2 + 2u (0) 



Z 



4((z-u(°)) 2 -4u(°)) 5 

^(o) n (o) 2 + s (z _ u (o) )w (o) n (o) + {z _ u (o) )v (o) u (o) + 3^(0)2 + 2v(0)v (o) 
H 5 

(( 2 _ n (0))2_ 4u (0))2 

10(* - V (°) W <?> + 10^°)\i 0)2 + 10 



+ 

which leads to 



4 24(^ 2 -,(0)^°) 2 ) 



+ 



24(4°) 2 -,;(0)nl 0)2 ) 2 



24(,f 2 -,(0)nf 2 ) 2 

_ ,W(^° )4 - 3ug^g + 2u< W> tjgj + A + 2(t,<°> ujf + ~ ^ tft) 

24a° )2 -,(0)n(°) 2 ) 

A further coefficient can be easily computed by taking j = 2 in (|66h . Thus we obtain 

u (3) = I u (2) _ _L„(0) (3) = o 

2 24 

4.2 The classical limit 

In the classical limit e — » the functions reduce to the first terms (u^°\ v^) of their 

semiclassical expansions and verify the equations of the dispersionless 1-Toda hierarchy 

d tj u = - d x (r j+1 - urj), d tj v = vd x rj, (70) 

where r,- are the coefficients of the Laurent expansion of R := 

R:= = n N9 a = 2^ fe — ' r o = 1 > 71 ) 

p-p v/(5 - u 2 - 4w ^ 



16 



and we have taken into account (see (|6ip ) that T = = (z — u) R/z. Here p := p^ and p := p^ 
are given by 



P{z) = \({Z-U) + y/( Z ~ u) 2 - 4v) 



Z — U h 

z 



P{z) = - {{z-u)- y/(z- u) 2 - 4vj = - 

According to (f64"j) . (u,v) are determined by 

dz U z 



+ 



2vri ^(z - u )2 - 4v 

dz z U z 



0, 



_ L, 27ri y/( z - u) 2 - 4v 
which can be also expressed as hodograph type equations 

Ejlii ((tj+cj) - (tj + Cj 



-2x, 



)Ki-i = 0, 



L 2 E"ii((*i + Ci)-fe+c j . 



r-j = x. 



4.3 The planar limit of the Hermitian matrix model 

From (|50p the one-point correlator W(z) is given by 

1 <9S r 



1 1 C7§ 

= 7 + M r (e,z, l,t) - e ^ — -^( e ,z, l,t), 



3>1 



so that by using the first equations of Q and ([5T|) one finds 

oo 

W{z) = J2(j + 1) (f i+1 - 6 ) (ay + - ^) 

3=0 ' 



i=i 



where 



H (z 



i-2 - 



■ T (zi~ 2 : 



e 



(72) 



(73) 



(74) 



(75) 



We are going to show that the solution of the dispersionless 1-Toda hierarchy determined by (|74p 
and (|49p describes the planar limit of the Hermitian matrix model in the one-cut case where the 
density of eigenvalues 

p(z) = M(z) ^( z -a)(z-b), 
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is supported on a single interval [a, b\. As it is known (see for instance |33j-|34j) these objects are 
related to the first term of the large ^-expansion of W in the form 

W<°> = ~\v z {z) + mp{z), V(z) :=Y,hz k . 

k>l 

According to (|6ip . in the classical limit T = = (z — u) R/z and then from (|75p it follows 

-j OO -j OO -. oo 



3=1 



3=1 



3=2 



with x = 1 in all x-dependent functions. Due to (l49j) the first hodograph equation (|74|) implies that 
the first term in the last equation vanishes. Therefore the expressions for the density of eigenvalues 
and the end-points of its support provided the above solution of the dispersionless 1-Toda hierarchy 
are 



p(z) :- 



1 



(- 



v; 



2m\ ( z - a){z - b) 



) y/( z -a)(z-b), 

J ffi 



(76) 



a := u — 2 y/v, b := u + 2 y/v, 



where x = 1 in all x-dependent functions. Moreover, from ([73]), they are determined by the 
equations 

dz V z 



7 2vri _ a )( z _ b) 
dz z V? 



(77) 



/ 7 2vri ^(z-a)(z-6) 



-2. 



They coincide with the equations for the planar limit contribution to the partition function of the 
hermitian model 1 30 1 ■ 



5 Critical points and the double scaling limit 

As we have seen the characterization of (u, v) as semiclassical expansions relies on the determination 
of smooth leading terms (u^ , v^), which are defined implicitly by the hodograph equations 
(|64j) . However, near critical points the functions (u^°\ v^) are multivalued and have singular 
x-derivatives. Thus the semiclassical expansions are not longer valid and a different procedure 
must be used. In this subsection we indicate how the so called double scaling limit method (see 
for instance [35]) can be formulated in our scheme. 
To simplify the discussion we set u = and 

t 2j -i = Cj = 0, j> 1; tj = cj = 0, j> 1, (78) 

so that the Lax operator is of the form 

C = A + vA- 1 , (79) 
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and we are only considering the Toda flows associated with the even times t<ij- After eliminating 
M in (133 j) . one sees that the generating function U := T[_i] satisfies the identity 

v (U + U[_i]) (U + U[i]) = z 2 (V 2 - 1). (80) 
This leads to expansions of the form 

U = ^%, U = ^6 2fc C/( fc ). (81) 

j>0 z k>0 

On the other hand, the system ([55]) reduces to 

oo 1 

^ U * = --T-- f dzV z U = x, (82) 
j=i 47 ™ J-y 

where V = Ylk>l ^2fe z 2k . Thus , the solution v is found from ([80]) and ([82]) . In particular, the 
leading term implicitly determined by the hodograph equation 

H(t even ,v^) = x, (83) 

where 

H(t even ,v) :=--?- <f dzV s U<® = --^ / dz 



Given a general m-th order critical point v c := v c (t even ) satisfying 



dH 

dv 



d m ~ L H 



dv 



m—l 



0: 



d m H 



dv r 



7^0, 



the method of the double scaling limit introduces a new small parameter e and a new variable x 
given by 

e:= e A, x = J?(u c ) + e™ x, (84) 
and generates solutions to ([80]) and ([82]) of the form 



1 + U = ^e fc C/W. (85) 

k>l k>0 

To prove it, we first observe that ed x = e 1//2 so that ([80]) can be rewritten as 

u E e "(7^v U ^" U + E ^ fc Uaiu) =(z 2 -4^)U 2 -z 2 , (86) 

n>l fc+i=2n;fc,/>l 



and by substituting the expansions (|85p in this identity and equating e-powers one can express 
each coefficient in the form 

u {n) = E /2 zv c ( 87 ) 

r=1 (z 2 -4v c ) 2 
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where the coefficients G n ,r are differential polynomials in tfW, 1 < k < n — r + 1 and their x- 
derivatives . In particular 

G„,i = 2«W 

and the first few fA 71 ) are 



i 3 ? 

(z 2 — 4v c )2 (z 2 — Av c )2 

~ (2) = 2^ 2 ) | (3nWVjnW) 
(z 2 - 4v c )^ (z 2 - 4v c )* 



2 ^ n (3) , c 2 z(l2u« (6n( 2 ) + al^))+12 9ln( 2 ) + a| 



^(3) _ f ^ » + 



(z 2 -4w c )2 6(z 2 -4v c )2 
2v c 3 z (I0(u( 1 )) 3 + 5(^n( 1 )) 2 + 10n( 1 )0|^ 1 ) +^uW) 



+ 



j 2 -4U,) 2 



Notice that = U^(v). 

By substituting (j8"31) - ([8"51) in ([82j) we get the system 



<p dz V z U (j) =0, j = 1, . . . , m - 1, 



I 4v™ y 



(88) 



Since v c is a m-th order critical point of (|83p we have that 

dz 2j+r = °i i = l,...,m- 1. 

7 (z 2 — Av c ) 2 

Hence, in view of (|87p . the first m — 1 equations in (|88|) are identically satisfied while the remaining 
ones become 



For n = m we get the equation which determines the the leading contribution uW in the double 
scaling limit 

G m , m (u {1) ) = K m x, K' 1 := v c m i — 2m+1 ■ (90) 



2 



For example 



2, 2(3u( 1 ) 2 + diu( 1 ) N ) =^ 2 5; 



m 



m = 3, 2 f 10 (u (1) ) 3 + 5 (<9 £ u«) 2 + 10u (1) dju« + = K 3 x 
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For n > m + 1 the equations of the system (|89p characterize the coefficients ti^ for k > 2. 

The differential equations ([90]) for vP 1 ' are essentially the stationary KdV equations [36] . Indeed, 
from (|8"6|) and taking into account flSTJ) one gets (Gj := G^j, G^ := %Gj, ...) 

2v c £ G ^ G j- E G.G.+W 1 ) £ GiGj-Vc Yl G i G 'j = °> 

i+j=m—l i+j=m i+j=m—l i+j—m—1 

which, up to trivial rescalings, coincides with the equation verified by the coefficients of the ex- 
pansion of the resolvent diagonal R of the Schodinger operator 

RR " _ 2 (z 2 -u)R 2 - -R' 2 + 2 z 2 = 0, R = l + Y?±. 

3>l 
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